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Abstract 



U(n) Yang-Mills theory on the fuzzy sphere Sfj is quantized using random matrix 
methods. The gauge theory is formulated as a matrix model for a single Hermitian 
matrix subject to a constraint, and a potential with two degenerate minima. This allows 
to reduce the path integral over the gauge fields to an integral over eigenvalues, which 
can be evaluated for large N. The partition function of U(n) Yang-Mills theory on the 
classical sphere is recovered in the large N limit, as a sum over instanton contributions. 
The monopole solutions are found explicitly. 
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1 Introduction 

Gauge theories provide the best known description of the fundamental forces in nature. 
At very short distances however, physics is not known, and it is plausible that spacetime 
is quantized below some scale. This idea has been contemplated for quite some time, 
and received a boost recently due to the discovery that string theory naturally leads to 
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noncommutative gauge theories under suitable conditions, as explained in 1]. Gauge 
theory on noncommutative spaces has been the subject of much research activity in 
recent years, see e.g. j2j for a review. 

There is one major problem with most models of noncommutative gauge theories: 
their quantization is very difficult. A direct quantization leads to difficulties related 
to the notorious UV/IR mixing |3], see also jl] for some recent developments. On the 
other hand, the use of the Seiberg-Witten map PIE], which allows a formulation in 
terms of commutative quantities, yields Lagrangians which become increasingly compli- 
cated at each order in the deformation parameter. This seems to rule out perturbative 
quantization 0. 

The motivation behind this paper is to try to develop new tools for the quantization 
of gauge theories, taking advantage of noncommutativity. The idea is to make use of one 
very fascinating feature of gauge theory on (some) noncommutative spaces: It is possible 
to formulate the gauge theories in terms of Lagrangians which have no derivatives. 
Rather, the dynamical variables are essentially matrices Bi, and the action is the trace 
of products of these matrices. Hence these gauge theories are matrix models. The gauge 
transformations have the form Bi — > U^BJJ for unitary matrices U. The kinetic term 
is generated upon a shift Bi = Xi + A4, and the Ai become the usual gauge fields in 
the commutative limit. This is very interesting for the quantization, because the path 
integral can now be defined simply as the integral over the matrices Bi, as in a random 
matrix model. A promising strategy is then to first do the quantization in terms of the 
Bi fields, and then go to the "classical" variables A4. This is a bit in the spirit of [Zj. 

In general of course, things are still complicated: The actions are multi-matrix models 
with nontrivial interactions, and the integration over the Bi is highly nontrivial. It must 
be so, since they describe a nontrivial quantum field theory. Moreover, the matrices are 
infinite-dimensional for most spaces (such as for M.$). This latter problem does not 
occur on the so-called fuzzy spaces, in particular the fuzzy sphere Sjj |Hj- This quantum 
space is characterized by a deformation parameter which measures the size of "Planck 
cells", and reduces to the classical sphere for N — > 00. Moreover the rotation invariance 
under SU (2) is maintained, hence Sfj seems particularly well suited to explore this idea. 

In this paper, we will show that for pure U(n) Yang-Mills theory on the fuzzy sphere, 
the quantization can be carried out completely by integrating over the matrices Bi. This 
is achieved by collecting the Bi into a single hermitian matrix, subject to a constraint. 
Of course, Yang-Mills theory on S 2 is a rather simple field theory with no propagating 
degrees of freedom; however it does have nontrivial monopole sectors, and its quantiza- 
tion is not entirely trivial. We will calculate the partition function for U(n) Yang-Mills 
theory on Sfj in the large N limit, and recover the known result [HI EI] for the par- 
tition function on the classical sphere. Corrections of order — could be calculated in 
principle, but we do not attempt this here. The main message is the applicability of 
completely new methods to noncommutative gauge theory, and hence to their commuta- 
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tive limit. Moreover, our result strongly suggests that the "commutative limit" of pure 
gauge theory on the fuzzy sphere is smooth, which is not obvious in view of the UV/IR 
mixing effects in noncommutative field theories jSj- This was also found recently on the 
quantum torus [llj, however with very different methods. 

Another important message is that in the approach developed here, the construction 
of gauge theory on (some) noncommutative spaces can be simpler than on a classical 
space. In particular, there is no need to introduce nontrivial fiber bundles, connections 
and other mathematical structures in our approach: the monopole sectors arise auto- 
matically in a very simple way, and reproduce the correct classical limit. We explicitly 
calculate the gauge fields for all monopole configurations. 

This paper is organized as follows. After a brief review of the fuzzy sphere and the 
most basic facts about matrix models, we present in Section 13.21 the particular poten- 
tial to be used in this paper. We then show that its minima define a fuzzy sphere, 
and the fluctuations become gauge fields on this fuzzy sphere after imposing a suitable 
constraint. In Section l3~Kl the monopole sectors of the U(l) case are identified, and the 
gauge field for the monopoles is calculated explicitly. We then generalize the construc- 
tion to the U(n) case in Section which amounts simply to taking larger matrices. 
Section El contains the calculation of the path integral, which is the main application of 
our construction. Finally we make some simple observations on symmetries, correlation 
functions, and show how a small modification leads to gauge theory on the g-deformed 
fuzzy sphere. The technical part of the path integral calculation is postponed to the 
appendix. In general, the focus is on explicit calculations, keeping the formal mathe- 
matics to a minimum. The hope is that noncommutative field theory in general and at 
least some of the techniques developed here will eventually be useful for physics. 



2 The basic fuzzy sphere 

We start by recalling the definition of fuzzy sphere [%|IT2*] in order to fix our conventions. 
The algebra of functions on the fuzzy sphere is the finite algebra Sjj generated by 
Hermitian operators Xi = (xi,X2,X3) satisfying the defining relations 



The noncommutativity parameter Ajy is of dimension length, and can be taken positive. 
The radius R is quantized in units of A N by 



This quantization can be easily understood. Indeed (JH) is simply the Lie algebra su(2), 
whose irreducible representation have dimension N. The Casimir of the iV-dimensional 
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representation is quantized, and related to R 2 by (J2J) and (JHJ). Thus the fuzzy sphere 
is characterized by its radius R and the "noncommutativity parameters" N or Ajv- 
The algebra of "functions" S 2 N is simply the algebra Mat(N) of N x N matrices. It 
is covariant under the adjoint action of SU(2), under which it decomposes into the 
irreducible representations with dimensions (1) © (3) © (5) © ... © (2N — 1). The integral 
of a function / e Sfj over the fuzzy sphere is given by 

r 2 j m = ^-nmi w 

where we have introduced J , the integral over the fuzzy sphere with unit radius. It 
agrees with the integral J dQ on S 2 in the large N limit. Invariance of the integral 
under the rotations SU (2) amounts to invariance of the trace under adjoint action. It 
is convenient to introduce the dimensionless coordinates 

Xi = Xi/A N (5) 

which satisfy 

•, jV 2 — 1 

£ k^j = *Afc, AiA 4 = — - — • (6) 

The Aj form a A^-dimensional representation of 577(2), which is given explicitly in Ap- 
pendix A for convenience. Noting that [A», xy] = iEijkX k , it follows that the rotation 
operators Jj act on functions / G Sj^ as 

Jif=[\,f]- (7) 
One can now write down actions for scalar fields, such as 

S = |i<f(A + ^)$ + |$ 4 (8) 

were $ is a Hermitian matrix, and A = ^ Jf is the Laplace operator. For gauge fields, 
the "correct" action is less obvious because the gauge fields have a priori 3 components 
(because there are 3 independent one-forms on Sjj, 0), and it is not obvious how to 
get rid of the normally unwanted 3rd component. Several slightly different approaches 
have been pursued in the literature jHlHSllIlllISllini • Alternatively, if one keeps the 3rd 
component which is essentially a scalar field, one finds actions which turn out to describe 
D2-branes on SU(2) jTHj. 

In this paper, we will develop a particularly simple formulation of gauge theory on 
Sjj, which makes a clear choice of the preferred actions even in the nonabelian U(n) 
case, and includes the topologically nontrivial sectors in an extremely simple way. The 
starting point is the following observation: we can combine the generators Aj which are 
N x N matrices into a single 2A^ x 2A^ matrix by 

C=\ + J2 X ^- (9) 
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We then observe the following property 



C 2 = (-) 2 , (10) 

which follows from (|139j) . Hence the eigenvalues of C are ±y. To get the multiplicities, 
note that AjO"* is an intertwiner of (2) <g> (N) = (N — 1) © (N + 1) (i.e. it is invariant 
under £77(2)), hence the multiplicities are N + 1 resp. N — 1. 

This simple observation leads to the idea that one should consider a matrix model 
for a hermitian matrix C = | + £?jcr l , and a potential which has ±| as degenerate 
minima. The fluctuations Bi = Aj + A, around the above solution should correspond 
to the gauge fields Aj, and the invariance under C — > U~ l CU for a unitary matrix U 
should correspond to gauge transformations (and other symmetries). Indeed this idea 
works. Before working it out, let us briefly recall some basic facts about matrix models. 



3 Matrix Models and the Fuzzy Sphere 
3.1 A brief review of single- matrix models 

We briefly recall some basic facts about matrix models which have found many applica- 
tions in physics. We refer the reader to [TTHI2T)] for excellent reviews and more references. 
Consider the matrix model of a single N x N hermitian matrix C with potential V(C). 
The partition function of the model is defined by 

Z = f dCe~™^ = I f[d Ci A 2 (c)e-S*^) (11) 
J •* i=i 

where Cj are the N eigenvalues of the Hermitian matrix C. Here 

He) = II( C * - <*) (12) 

i<j 

is the Vandermonde determinant, which is the Jacobian of the transformation dC = 
Y\ dcidU A 2 (c) , and the integral over the unitary matrices U is trivial. In the literature 
on matrix models, the potential V is usually chosen to be of the form 

V(C) = - 2 v(C), v {C)=Y,9kC k (13) 

" k>2 

where #2 = 1 and the couplings are kept fixed in the large iV limit. 

The reason why these matrix models are so useful is that the models really only 
depend on the N eigenvalues q, while the matrices have iV 2 degrees of freedom. This 
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lead to the development of powerful methods (e.g. steepest descent method, orthogonal 
polynomials, etc [2HE2]) wrucn can be used to analyze the models, and basically solve 
them explicitly in the large N limit. For example, the saddle-point equation is given by 



1 l( , 2 



g 2 ' v ~* y N ^ a-Cj 

The sum in the r.h.s. is due to the Vandermonde determinant in the measure and rep- 
resents a repulsive potential among the eigenvalues. The iV-dependence of the potential 
(jlrS)) ensures that the repulsive effect is well balanced in the large N limit. Due to 
this repulsive force, the eigenvalues spread evenly around the classical solution of the 
equation of motion q = 0. The distribution of the eigenvalues 

P(c):=^£<Kc-q) (15) 

i 

becomes continuous in the large N limit and can be solved easily from the equations 

1 ./^ o f ^P(C) 
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v\c) = 2 J- dc'-^jj, j dcp(c) = 1. (16) 



There is much more to be said about these matrix models. In particular, the distribution 
of eigenvalues (and correlation functions) can be derived using e.g. the method of 
orthogonal polynomials, without relying on the saddle-point approximation. However 
we will not need these techniques due to the simplicity of the model considered here. 



3.2 A matrix model with degenerate minima 

In this paper, we will show that the fuzzy sphere arises as a vacuum solution of another 
Hermitian matrix model, given by a potential with a different scaling dependence in N. 
Consider the matrix model with action 

S = ^ = ^ (< C2 - 4 )2)2 ) - ^N Tl (£ " ^ + (17) 

where C is a 2N x 2N Hermitian matrix, and g 2 > is kept fixed independent of N. 
The shape of the potential is sketched in Figure [T] The following features distinguish it 
from the matrix models considered before: 



• The coefficient of the quadratic term is negative. This implies that the distribution 
of eigenvalues are peaked at the minima ±y of the polynomial V(C), rather than 
around the origin. 
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Figure 1: V(C) 



• The specific form and the particular N dependence of V(C) is chosen such that the 
accumulation of eigenvalues at its minima will lead to the emergence of a (fuzzy) 
sphere Sfj. The fluctuations will describe a gauge theory on Sfj, as we will show 
below. Apart from the N dependence, the relative coefficients between the C 2 and 
C 4 terms in V can be adjusted to any (negative) number by a rescaling of C and 
g. In general, different spaces may be generated for different potentials V(C), and 
the properties of this space such as symmetries are related to the details of this 
eigenvalue distribution. 

• If we expand V(C) around one of its minima (consider y to be specific, setting 
C = y + /i + ), then 

+^ + ) = ^(/4 + |/4 + ^/4) (is) 

Now the coefficient of the quadratic term is positive, and comparing with the usual 
matrix models (fT3|) we expect that the higher-order terms will become irrelevant 
for large N, leading to a simple Gaussian distribution near each minimum. 

Let us turn to the eigenvalue distribution in the large N limit. The stationary points 
of our potential are given by the eigenvalues 

N , \ 

c=0, ±y (19) 

Since the action for c = is of order N 3 , it is quite obvious that this stationary point 
does not contribute for large N. Hence consider the eigenvalues distribution near the 
minima c = ±y in the large N limit. We introduce q = y to facilitate the analysis; 
the other case q = — y + //j is similar. Performing a similar analysis as 21\, the saddle 
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point equation becomes 



1 , 6 9 4 . 1 2 

_(^ + - rf + _ ft ») = -g_. (20) 



In the large N limit this becomes again 
which gives the standard distribution for Gaussian matrices 



-// = / d//'^;, / d/ip(/i) = 1 (21) 



a = , p(ji) = — vV - V?- (22) 

g-K 

In particular, the distribution is only nonzero for < V%g, and we find a finite spread 
of eigenvalues. The reason is the explicit N in front of the quadratic term in the action. 
This means that only eigenvalues near ±y will contribute, 

6/n < V2g. (23) 

Hence near each minimum ±y, the action can be replaced by a Gaussian ^-^(/-t^) 2 
for large N, since the higher-order terms are suppressed by -L. 



3.3 Matrix model vacua and emergence of the fuzzy sphere 

Let us resume the analysis of the particular model given by the polynomial ([17)1 . The 
solutions to the classical equation of motion 

C{-N 2 /A + C 2 ) = (24) 

are characterized by the multiplicities n + ,n_,no of the eigenvalues ±y resp. 0, which 
satisfy n + + n_ + no = 2iV. We can assume that no = as discussed above, since each 
zero eigenvalue gives a contribution to the action and is highly suppressed. Then 
the saddle points are characterized by the trace Tr(C), given by 

Tr(C) = N/2 (n+ -n_). (25) 

Consider the vacuum with n + — n_ = 2, i.e. n + = iV + 1, n_ = N — 1 . Using the U(2N) 
invariance, one can put the vacuum in the form 

C = l + A»er' (26) 

where Aj are precisely the N x N Hermitian matrices in (JSJ, (JHJ) which describe the fuzzy 
sphere, and a 1 , i — 1, 2, 3 are the Pauli matrices. In other words, we can find a unitary 
matrix U such that 

diag(^, ....,—) = U{\ + X i a i )U' 1 (27) 
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provided n + — n_ = 2, as explained in Section El The equations of motion (J%4*|) then 
take the form 

£ k^'^j = i^ki ^i^ 1 — — ^ — • (28) 

These are precisely the defining relation © of the fuzzy sphere Sfj in terms of the 
dimensionless coordinates. The radius has been set to unit here, since it can easily be 
reintroduced. 



3.4 Matrix fluctuations and gauge fields 

Now consider a general 2N x 2N Hermitian matrix C, 

C = C a a a = (^ + p)cr + B l a i (29) 

where a = 1 . Plugging this into (|17jl. we obtain 

S = TrV(C) = Tr((5^ - A.A') 2 + (B t + ie ijk B 3 ; B k )(B i + ie irs B r B s ) 

+D i pD i p + N 2 p 2 + 2p 3 + p 4 

+6(5^ - A;A>(p + 1) + Aipe ijk (B' l B ] B k - A i A J 'A*)Y(30) 

where 

Ap:= (31) 

This starts to look like a field theoretic action on the fuzzy sphere Sfj. Its interpretation 
is however obscured by the presence of p. Comparing with , we shall therefore impose 
the constraint 

P = 0, i.e. C = ± (32) 

This implies 

Tr(C) = N. (33) 

Then the above action becomes 

9 / /V 2 — 1 \ 
5 = -^Tr^BiB* — ) 2 + (B, + ie ijk B*B k ){B i + ie irs B T B s )). (34) 

This is one possible action for a gauge theory on the fuzzy sphere, cp. |8lll3lll4pi5|ll6|ll7| . 
The constraint breaks the original SU(2N) symmetry down to a smaller subgroup, which 
contains a SU(N) gauge symmetry acting as 

Bi -> U^BiU. (35) 

We will see that this corresponds to the usual U(l) local gauge symmetry in the classical 
limit. Of course, breaking the full SU(2N) symmetry by the constraint (}32|) is somewhat 
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against the spirit of the matrix model, in particular the integration in the partition 
function cannot be carried out as easily as in (jllj) any more. This is to be expected, 
since f!34|) is a matrix model with 3 interacting matrices. Nevertheless, this approach 
will allow us to carry out the path integral with some more effort. 

It is easy to understand the dominant configurations for the action (|3*4*j) . The term 
(Bi + iEijkBi B k ) 2 implies that the Bi approximately generate a representation of su(2), 
and the other term implies that BiB % « (N 2 — l)/4, which corresponds to the Casimir 
of the approximate su(2) representation. Hence one could interpret ()34)) as a theory of 
"fluctuating representations" of su(2), and the dominant configurations will be approx- 
imately iV-dimensional irreps of su(2). This is an important difference to other possible 
actions without the term (B i B' t — JV 4 ~ 1 ) 2 , such as in jSH]: there, reducible "block"- 
solutions with blocks of arbitrary size are allowed, while in ()34)1 they are suppressed. As 
we will see, this is crucial for the physical interpretation, and only (jHIJ) reduces to an 
ordinary Yang-Mills theory on S 2 in the large N limit. 

Consider next the equations of motion, 
iV 2 -l 

[B\ BjB 3 — ]+ + (B + ieBB)i + ie l]k [B h (B + ieBB) k ] = 0. (36) 

The "vacuum" solution is 

Bi = Aj, (37) 

up to gauge transformation. In fact then S = 0, and this is the unique solution with 
S = up to SU (N) gauge invariance, because both B + ieBB = and B^B 1 - 
must hold. This means that Bi is a representation of su(2) with fixed Casimir 
now expand 

B\ Aj -f- ^4j, 

then 

B t B { - A i A i = \iA { + AiX* + AiA* 

and 

B< + = \e*>F u , F U := UK, A] - A k ] + i[A k ,M + s Um A- 

Notice that the kinetic terms in Fm arise automatically due to the shift (|38|). The SU(N) 
gauge symmetry acts on Ai as 

At -> U^AfU + U- l [\, U] (41) 

which for U = exp(ih(x)) and N — > oo becomes the usual (abelian) gauge transformation 
for a gauge field. 

However, the gauge field Ai has 3 components, which does not seem to match with 
the degrees of freedom in a 2-dimensional gauge theory. To understand this, we should 
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N 2 -l 
4 

If we 

(38) 

(39) 
(40) 



Then 



translate the action into conventional field theory language, which can be done for large 
N. One can then decompose the field Ai into a tangential component A\ and a radial 
component ip as follows 1 : 

4 2 

Ai = N2 _ Xiip + Ai « — Xi ip + Ai (42) 

where ip and Ai are defined such that 

XAi = 0, (43) 
p := XAi. (44) 

- A,A* = 2<^ - [A*, ^] + Aj + 1 T(^, A). (45) 

Here T(</>, A) stands for functions of (p and Ai which are suppressed by 4. Similarly, all 
terms involving ip in the "field strength" F k i (|40|) are suppressed 2 by Therefore the 
only term involving ip which contributes for large N is the square of ()45j) . We can now 
simply integrate out <p (i.e. consider it an auxiliary variable), replacing it by 

<P = l(lWA t ]-A*) (46) 

for large N, which is smooth (i.e. high angular momenta are suppressed, assuming the 
Ai are smooth). Hence all terms in F k i containing ip can be omitted for large N, being 
suppressed by jr. Another way of arriving at this conclusion is to use the field <fi := jyip, 
which has a large mass of order N. At any rate, we can now write 

F kl = i[A fc , Aft - i[X h A k ] + i[A k , A t ] + e klm A m (47) 

for large N, involving the tangential gauge field only. As a consequence of (|4*3^) and ©, 
F k i is "tangential" 

x k F kl = o(l/N) (48) 
for large N, and it becomes the field strength of an abelian gauge theory on S* 2 . 

To summarize, we found that the radial fluctuations ip decouple in the large N limit, 
and (jH3J) reduces to a U(l) Yang-Mills theory on a unit sphere with action 

S = \ J F mn F mn . (49) 

Here 



kl 



iJ k A t - UiA k + e klm A m (50) 



this decomposition as defined here is gauge- invariant only in the large N limit. From that point of 
view if := BiB 1 — AjA 1 would be a nicer radial field, however this leads to a nontrivial Jacobian in the 
path integral, which makes the decoupling argument below more subtle. 

2 assuming that <p and Ai is "smooth" , so that [<p, Xi] is finite. This is justified by (j46(l and the kinetic 
terms in the action. 
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is the field strength for the U(l) gauge potential A±. The fields are tangential in the 
sense [21] 

iMi = 0, x% = 0. (51) 

This is a description of the 2d gauge theory in terms of a 3-component gauge field A4 
subject to the tangential constraint (|51)l. This formulation is manifestly invariant under 
£0(3) rotations. To put it in a more familiar form, let us assume that we are sitting on 
the north pole. Then only Ai,A 2 and F 12 survive the constraints, and 

zJi = -d 2 , U 2 = d x . (52) 

Our gauge theory can then be identified with a gauge theory with only tangential gauge 
fields A\ d \ i = 1, 2, whose field strength takes the usual form 

F[? = d l A 2 cl) - d 2 A[ d) (53) 

if we identify 

A { f ] = -A 2 , 4 C ° = ^i- (54) 
In coordinate independent form, this is 

A id) —fx A (55) 

where f is the radial unit vector. This identification will be useful in the next subsection. 

Since the volume of the gauge group is finite here, we do not have to fix the SU(N) 
gauge using e.g. the Faddeev-Popov method. Instead we can keep the integral over all 
configurations. Indeed, working with Bi or even C seems much easier, and makes all 
the symmetries manifest. The beauty of our formulation is that it allows to apply the 
powerful methods of random matrix theory, after suitable modifications. One can hope 
that this will lead to new methods for studying gauge theories. 



Alternative version of constraint. It is important to realize that the matrix model 
TrV(C) describes a Yang-Mills theory only if we impose the constraint ()32j) . so that the 
last terms in ()30|) vanishes. Without that constraint, Aipe i jk{B l B : ' B k — XV \ k ) contains 
a term oc NFp, which after integrating out p cancels the YM term TtFF. However, 
there is another possibility, namely to consider 

3' = Tr(V(C)) - NTr(C - i) 2 . (56) 

The last terms of course implies the Co = | constraint in the large N limit. This leads 
to an additional term —TrNp 2 in the action (jHUj) . which allows to integrate out p for 
large N leaving a rescaled YM term oc TrFF. The path integral can then be carried 
out in the same way as we will do below. Basically this seems to be a matter of taste, 
and we chose to impose p = directly. 
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3.5 Monopole sectors 



If we claim to have a fuzzy version of Yang-Mills theory on the 2-sphere, we should be 
able to recover the monopole sectors as well. They are not hard to find here: as discussed 
above, the dominant contributions for the action (|34|) are "approximate" representations 
of su(2) with Casimir B,iB % « N 7 1 ■ This suggests to consider irreps of dimension M 
slightly different from N . They indeed turn out to describe monopoles 3 . 

Hence consider the same action p7jl . but for matrices of different size. Let 

= \ + B\ M) a* (57) 

be a 2M x 2M matrix with 

M = N — m, m G Z. (58) 

This implies in particular 

Tr(C) = M, (59) 

which again picks out the sector n + — n_ = 2. One can easily see that the equation of 
motion (|3T)j) resp. Tr(V'(C)5C) = has solutions of the form = a m X[ M ^ if 

3{a m - l) 2 + a 2 m M 2 - N 2 = 0, (60) 

which gives 

Tfl 

a m = 1 + — for > m. (61) 



Hence we found the new solution 



C^ = \ + a m ^a\ (62) 



Then 
and 



(F (M) )i = fej'W")* + (S (M) ), - y x t (63) 

N 2 - 1 

5 ( a/) ■ B {M) > (64) 

for N — > oo. In particular, the field strength is tangential in the sense ()48j1 . with 
|F| = y. This is just like the field strength of a monopole of charge m, suggesting that 
m is the monopole charge. The action is 

S(CM) = g (65) 



3 The basic idea that fuzzy spheres of different size correspond to monopoles was also proposed 
in [23, without calculating the gauge field. 
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for large N. We will show that this interpretation as a monopole is correct by writing 
as an excitation over the fuzzy sphere solution Bi = \\ . The corresponding 
gauge field will take the usual form of Dirac monopole of charge m in the large N limit. 
To see this, write this solution block matrix 

1 , f a^\\ M) ^ 0\ 1 , ,(N) 



C =2 + { oJ = 2 +A ' rV + ^' (66) 
where Ai is interpreted as gauge field. Hence 

A, = a m \^-\f\ (67) 

Using the representation ()136|) for the A, = , one obtains the following non- vanishing 
matrix elements 

= \ ~\"t (1 - Q ™ (1 " ) I N ~ m ' (68) 

I -{M)kh, k > N — m, 

(A + ) kM1 = i i 1 ~ «™>A=^ ) > 1 < ^ < ^ - m, (6Q) 

[ -(A+)fc,fc+i, k> N -m 

and A_ is obtained by transposition. This can be translated into functions on the fuzzy 
sphere via (0), which takes the form 

Xi = |a< (70) 

in the large N limit. We also note that k is related to the "height" x 3 through 

N + 1 

k = - (AsV (71) 

The quantities in ([68 )1 - ([69 )1 have a smooth limit for large N, except at the finite set of 
"points" N — m+1 < k < N (located at the south pole) where Ai develops a singularity. 
This singularity corresponds to the Dirac string. In the patch covered by 1 < k < N— m, 
which represents the sphere without the south pole, we obtain 

A 3 = --(l-x 3 ), (72) 

mx- mx- 

A- = 7 r H , (74) 

2(1+ s 3 ) 2 K J 

in the large N limit (recall that R = 1 throughout). It is easy to check that Ai satisfy 
the constraint (|5T]) . This looks almost but not quite right; however, recall that we must 
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use the identification (|55j) to find the corresponding classical gauge field Af\ This 
comes indeed out as 

A^ = rxA = ^-^— ( X | (75) 



2 l + x 3 







or = y (ydx — xdy), which is precisely the (tangential) gauge field of a Dirac 
monopole of charge m on the sphere. The field strength was already calculated in 
and is constant with the correct quantization 

/TYl 
F iX l = 4tt- (76) 

Notice that F is constant in spite of (or rather because of) the "non-classical" term 
[Ai, Aj] in the definition of F. The same calculation applies for m < 0, hence we get 
both negative and positive monopole charge as it should be. The singularity at the south 
pole can of course be moved around using suitable SU(M) gauge transformations. At 
finite N resp. M, this configuration should therefore be interpreted as a fuzzy monopole. 

This point of view considering the monopole sectors as matrices of different size 
is quite compatible with the treatment of nontrivial topological sectors in j2H|, where 
sections in nontrivial bundles are represented by iV x M matrices. Clearly our gauge 
fields Bi of the appropriate size can act on these from the left resp. right, and one can 
define covariant derivatives in this way. This will be elaborated elsewhere. 



More careful embedding of monopole sectors. We should address a somewhat 
unsatisfactory aspect of the above treatment of the monopole sectors: Formally we have 
been considering distinct matrices for different M, but more properly they should 
all be considered as block-matrices embedded 

... ^ Mat(N -1)^ Mat(N) ^ Mat(N + 1) ^ ... (77) 

as in (JBTIj) . cp. also j2ZJ. Then there is a small problem with (JBT)|) : the | in the lower- 
right block in ()66|) must be there in order to satisfy the constraint p = 0. However the 
eigenvalues of this small block are far from ±y. Therefore this type of block-matrix 
configuration would strictly speaking be highly suppressed by the action, because the 
Dirac-string contributes o(N 3 ) to the action. One way to cure this problem is to replace 
the action by 

S' = Tr(^i^V(C)), (78) 

which now has ±y, | as degenerate minima, and all fluctuations are Gaussian as in 
(jl8j) . Here C G Mat(M) should be a matrix of fixed size M which is large enough to 
accommodate all relevant solutions, i.e. IfcMc 2N. Now the block-matrix (}66|) is 
really a solution of the equation of motion, with action ^2. In terms of the Bi fields, 
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we note that (C — |) 2 = B i B i + ie lkl B k Bia k , so that the action differs from by- 
terms of the form Tr((-^- + -^Xi(7 l )y(C)) which are suppressed for large N. The crucial 
difference is that the Dirac-string now has vanishing action. Therefore all monopole 
configurations can be obtained as distinct solutions in one single configuration space 
Mat(M). This is conceptually very appealing, because it shows that the nontrivial 
topological sectors arise here automatically as different solutions for the same action. 
This is even simpler than in the classical case: there is no need to introduce nontrivial 
principal bundles, they just come out. However the calculation of the partition function 
below would be somewhat more complicated for the action (|78|). Since the path integral 
is the main focus of this paper, we shall not pursue this point of view here, and consider 
the monopole configuration as truly "distinct" sectors for simplicity, as classically. 



4 Nonabelian case: U(ri) Yang-Mills theory 

Now consider the same matrix model S = TrV(C) as in ([17)1 . but for larger matrices of 
size 2M x 2M with 

M = nN — m, Tr(C) = M. (79) 

The last constraint implies that the multiplicities of the (dominant) eigenvalue distri- 
butions of C are now n + — n_ = 2n. We will see that this leads to a non-abelian U (n) 
Yang-Mills theory. 

First we should find the ground state. For M = nN, the absolute minima of the 
action S = TrV are now given by any matrix C with n + = M + n eigenvalues +y and 
n_ = M — n eigenvalues — ^. In a suitable basis, C takes the form 

C = {\ + AS*V) l„xn, (80) 

which is a block matrix consisting of k blocks of the solutions (| + X^a 1 ) of Section 
13.31 The action is then zero, and clearly all other saddle points have a positive action. 
In general, we can write again any 2M x 2M matrix C in the form 

C=(± + p) + B i o i (81) 

where the Bi and p are now M x M matrices. In order to obtain a Yang-Mills gauge 
theory, we shall impose again the constraint 

P = 0, (82) 

so that the action (fTTj) reduces to 

2 / N 2 — 1 \ 
S = ^ Tr (( 5 ^ J—)' + ^ + i£ijkB i B k )(B i + ie irs B r B s )j, (83) 
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which has the same form as ()34j) but for different size M of the matrices. It is invariant 
under the gauge transformations Bi — > U^BiU for U G U(M). To understand its 
meaning, we write the fluctuations of B resp. C in the form 

B t = B h J a = xf ) t° + At (84) 

where Bi and carry a index, 

A = A ifl t° + A i:a t a . (85) 

Here t a denote the Gell-Mann matrices of su(n), which satisfy 

tJb = ^9ab + ^(d ab C + ifab C )t c (86) 

and to = 1 is the nxn unit matrix. The rest of the analysis of Section l3~4l goes essentially 
through. In particular, we can split the gauge fields again into tangential and radial 
components 

4A- 

We suppress here the u(n) labels of (p and Aj, which are defined by 

(p := X i A i , \ l Ai = 0. (88) 

Then (145)1 implies as before that all components of ip decouple and can be integrated 
out. It remains 

B t = Xi + Ai = Xi + A ifi t° + A iA t a (89) 

involving only the tangential components of Ai, and in the large TV limit we obtain a 
theory with action 

S = TrV = ^J F mn F mn (90) 

where again 

B' + is-B k B, = \e*>F tl , F k , = i[X k , Ai] - ,[*,, A k] + i[A t , A,] + wl- (91) 
is tangential 

x k F kl = o(l/N). (92) 

Spelling out the u(n) structure explicitly and omitting terms which vanish for large N, 
this action becomes 

S = TxV = j 2 J (F mn , F mn >° + F mn , a F m "< a ), (93) 

where 

Fki,o = iJkAifi — iJiA k0 + Eki m A mj0 , 

F k l,a = iJkAl^ a - UlA k)a + lA^A^f^ + e k l m A m)a . (94) 



This is the action of a U(n) Yang-Mills theory on the sphere. Recall that the only 
difference to the abelian case in Section 13.41 is the size M « nN of the matrices. 
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Saddle points. The remaining saddle-points can now be found from the equation of 
motion (J3H|) as in the previous section. Clearly any (reducible, in general) representation 
of su(2) with a suitable normalization as in (|62j) will give a solution. Therefore the 
(dominant) saddle-points are given by the block-matrices 



/ 






where each block has the form (|fi2|) with size 








\ 



(95) 



Mi = N — rrii, rrii € Z, 



such that mi + ... + m n = m. Their action is 

g((j(mi,...,m n )\ _ 



— y 

2,f- 4- 



mf 



(96) 



(97) 



for large N, hence configurations with large \rrii\ are suppressed. 

We can now write the saddle-points ([9*3|l as fluctuations around the ground state, as 
in Section 1331 After arranging the blocks appropriately (by a gauge transformation), 
they take the form 

/ m x Ai ... \ 
m 2 Ai ... 



.4; 



V 



o 



o 



Tfl n Ai 



(98) 



/ 



(for large N) with mi + ... +m n = m, where A{ is the basic abelian monopole field found 
in Section 13.51 Hence the sectors with rrii 7^ correspond to nontrivial U(n) gauge 
field configurations, which are precisely the "instantons" 4 of the U(n) YM theory found 



in 



5 The path integral 



The quantization of gauge theory on the 2-sphere has been studied extensively, us- 
ing a variety of methods including lattice formulations and a generalization of the 
Duistermaat-Heckmann localization theorem, see e.g. 01101 12E HIE USB EH In 
particular, the partition function and correlation functions of Wilson loops have been 

4 we refer to any critical point of the YM action as instanton, as is customary in the related literature. 
In general they are unstable. 
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calculated. It is therefore natural to ask whether such calculations can also be done on 
the fuzzy sphere. Some of the known methods, in particular the localization theorem, 
might well be applicable as shown in ^T] for the case of a torus. However this method 
is rather indirect, and we want to calculate the path integral directly taking advantage 
of the above formulation as matrix model. 

One of the nice features of the fuzzy sphere is the fact that all path integrals are 
finite, simply because there are only finitely many degrees of freedom. However, this 
does not necessarily make them easy to evaluate: e.g. for scalar fields, one is forced to 
resort to perturbation theory (see e.g. [51]). which is even more complicated than in the 
classical case. 

The main advantage of our matrix formulation of gauge theory is that it allows to 
explicitly carry out the path integral. This provides a truly new approach to gauge 
theory, since this formulation is possible only in the noncommutative case. Without the 
constraint p = 0, the integration would even be "trivial" as in Section ITT] but this does 
not describe a YM theory on the sphere. We will now show how this constraint can be 
handled using the known matrix model technology, and calculate the partition function 
directly by integrating over the gauge fields for large N. We will recover the known 
result |H1E3] for the partition function of U(ri) YM theory on the sphere for N — > oo. 
While the explicit calculation in Appendix B may seem a bit involved for our present 
application, one can hope that the idea will be useful in less "trivial" cases as well. 

We want to quantize the gauge theory with action by integrating over the 
M x M matrices Bi. This will be done by integrating over the 2M x 2M matrices C 
in the action (j!7j) . imposing the constraint (|32j) . We will not attempt here to calculate 
the full generating functional for the gauge field, only the partition function 



dBi exp(-S(B))) 
= j dC 5(C - ~) exp(-TrV(C)) 

= J c/AjA 2 (Aj) exp(-W(A)) J dU8{{IJ- l MJ) G - ^) (99) 

where dll is the integral over 2M x 2M unitary matrices, and C = U~ 1 AU. Here 
5{Cq — |) is a product over M 2 delta functions, which can be calculated as follows: 
define 

J=( K n l)=K*° (100) 



K 

where K is a M x M matrix. Then 



5{{U- l CU) Q - ~) = j dKeM^{U-\C-^)UJ)). (101) 



19 



By gauge invariance, the r.h.s. depends on the eigenvalues Aj of C only. Hence 
Z = J dK J dA i A 2 (A i )exp(-W(A)) / dU exp(zTr(£/ _1 A£/J — i J)) 

dK Z[J] e"2 TrJ (102) 



where 



Z[J] ■ = J dC exp(-TW(C) +iTr(CJ)) (103) 
depends only on the eigenvalues Jj of J. Diagonalizing K = V~ l kV , we get 



dhA 2 (k) J c/A i A'(A i )exp(-TrF(A)) J dU ex-p(iTr(U~ l (A- -)UJ)) (104) 

where J dV was absorbed in J dU. The main step is now to carry out the integral over 
J dU, which can be done using the Itzykson-Zuber-Harish-Chandra formula [SEIEE], 

/det(e iAiJj ) 
dUexpiiTriU^CUJ)) = const A ^^J^ • ( 105 ) 

This depends only on the eigenvalues of J and C, with Vandermonde-determinants 
A(Aj) and A(Jj). Note that the Vandermonde-determinants are totally antisymmetric, 
and so is det(e* AiJj ). Therefore this expression is manifestly symmetric in both A and 
J. 

In this step we have reduced the number of integrals from M 2 to 2M. This means 
basically that the integral over fields on Sjy is reduced to the integral over functions in 
one variable. This is a huge step, just like in the usual matrix models. The constraint 
however forces us to evaluate in addition the integral over hi, which is quite complicated 
due to the rapid oscillations in det(e tAtJ: >); recall that Aj « ±y. Nevertheless, it is 
shown in Appendix B how the integrals over Aj and ki can be evaluated for large N, 
with the result 

/oo 2 
dK X ...dK n A 2 (k) eX p(-|-^K 2 ). (106) 

for matrices of size M = nN — m (omitting overall constants). This form of Z was found 
in (HI] for a U(n) Yang-Mills theory on the ordinary 2-sphere, apart from the constraint 
^rrii = m which will be removed soon. It can be rewritten in the "localized" form as a 
weighted sum of saddle-point contributions, as advocated by Witten |2H]- This can be 
seen as follows: 



m\-\-...+m n =m 



20 



m\-\-...+m n =m 



"dm; 2g 

m\ + ...+m n =m 

E p ( m ^) ex p(-iE m ^)- ( lor ) 

mi+...+m«=m 

Here P(rrii,g) is a totally symmetric polynomial in the mj. The last exponential is 
precisely the action (J9?") for the saddle-point (mi, m n ) as discussed in Section[~J which 
is weighted by the polynomial P(rrii,g) (e.g. for n = 2, one finds P(jrii,g) = (mi — 
7B2) 2 — 2g 2 ). This shows that the "localization" |28j also holds in the noncommutative 
case for gauge theory on the fuzzy sphere, at least in the large N limit. However we did 
not use any localization theorem here, it comes out by an explicit computation of the 
purely bosonic path integral, without having to introduce auxiliary fermionic fields as 
in j2Hj- Note in particular that we did not do any gauge- fixing, which is not necessary 
here because the volume of the gauge group U (M) is finite. 

To include all monopole configurations, we should sum over matrices of different 
sizes M = nN — m as explained in Section 01 keeping V(C) constant. Then m is the 
U(l) monopole charge. Hence the full partition function is obtained by summing 5 over 
all Z m , 

,2 



/* 2 

z = Y,z m = J d * a2 ( k ) elKimi ex P(~y E < 



(108) 



m mi % 

One can now perform a Poisson-resummation as in [""I" . 



Y,f(m) = J2f(^p) (109) 

m p 

where f(p) = J ^f(x)e~ ipx . This gives 

Z= A 2 (p) exp(-27rVE^)- ( 110 ) 

pi,...,p n ez i 

This is the partition function of a U(n) Yang-Mills theory on the ordinary 2-sphere. As 
shown in [31], this is equivalent to the form found in piTHj 

Z = J2 ( d *) 2 exp(-47r 2 ^ 2 C 2i? ), (111) 

R 

where the sum is over all representations of U(n) and is the dimension of the represen- 
tation and C2R the quadratic casimir. Hence the limit iV — > 00 of the partition function 

5 the relative weights of Z m for different m is strictly speaking not determined here. However, it 
could be calculated using the embedding (|77|l as explained in Section 13*51 
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for U(n) YM on the fuzzy sphere is well-defined, and reproduces the result for YM on 
the classical sphere. This strongly suggests that the same holds for the full YM theory 
on the fuzzy sphere, and that there is nothing like UV/IR mixing for pure gauge theory 
on S%. This is unlike the case of a scalar field, which exhibits a "non-commutative 
anomaly" [33] related to UV/IR mixing. 



6 Remarks on symmetries and correlation functions. 

Let us try to understand the symmetries of our model in more detail. Recall that 
2-dimensional Yang-Mills theory is invariant under area-preserving diffeomorphisms 
(APD's), because the field strength F = Fijdx l dx 3 = fuj has only one component, 
and J F * F = J f 2 u is invariant under APD's. Here uo denotes the volume form. 

It is easy to understand the quantization of APD's on the fuzzy sphere 6 . The fuzzy 
sphere arises as quantization of the Poisson structure {xi,Xj} = EijkXk, which corre- 
sponds to the (canonical) symplectic form u = x l dx 3 dx k Eijk on S 2 which coincides with 
the volume form. Hence any function f(x) on the sphere defines a Hamiltonian flow, 
which preserves u. It is therefore an area-preserving diffeomorphism. Explicitly, the 
vector field Xf generating the APD with "hamiltonian" / is determined by 

£x f g = {g,f} (112) 

where {, } are the Poisson brackets. Hence the derivation 

Cx t (xi) = {xi, f} (113) 

is an infinitesimal APD determined by / acting on the coordinate function X{. After 
quantization, this is replaced by the commutator 

S f (xi) = [xi,f\ (114) 

which can therefore be interpreted as quantized infinitesimal APD on the fuzzy sphere. 
In this way, the APD's on the fuzzy sphere can be identified with SU(N). 

Now consider "abelian" gauge theory on Sfj as discussed in Section El It is well 
known (see e.g. [2B]) that gauge transformations on noncommutative spaces are closely 
related to certain diffeomorphism groups. Indeed according to the above discussion, the 
gauge transformations 

Bi -> U- l BiU (115) 

with U = exp(if(x)) G SU(N) can be interpreted as APD acting on the Bi, viewed as 
3 scalar fields. This statement has no classical analog, and has nothing to do with the 
action of APD's on classical gauge fields A4. The induced transformation 

Fij -> U- l F i3 U (116) 

6 I want to thank S. Rajeev for explaining this to me 
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of the field strength becomes trivial in the commutative limit, since then Fij commutes 
with all functions. For finite N however, 1)116)1 is nontrivial even in the "abelian" case 
at hand, and can be interpreted as action of the APD determined by U = e 1 ^^ on 
each component of F^ (which is again different from the classical action of APD's on 
a tensor!). Hence gauge transformations are necessarily "non-local", and cannot be 
disentangled from general coordinate transformations on S 1 ^. This phenomenon is quite 
common on noncommutative spaces: e.g. on the canonical quantum plane Kg, the 
translations are inner and hence part of the gauge transformations of F^. One could 
interpret this from a fiber-bundle point of view by saying that base space and fiber 
become unified in some sense. 

Apart from this SU(N) group of gauge transformations (or APD's), the constraint 
Co = | is also preserved by the SU(2) group exp(iaia l ) which acts on the indices of 
the Bi only. The overall symmetry group is therefore SU(N) x SU(2). Note that the 
"physical rotations" exp(i«j(A l + are combinations of these SU{2) rotations and 
gauge transformations. 

Let us extract some information about the correlation functions using these symme- 
tries, without trying to calculate them explicitly. We only consider the abelian case for 
simplicity. The correlation functions are defined by 

(C x ...C m ) = ^ J dC C x ...C m 6(C - l -) exp(-TrV(C)) (117) 

where the indices 1, m indicate 2N x 2N matrix labels, or in terms of the components 



((B l ) 1 ...(B,) m ) = 




dB (SOi-^Om exp(-TrV(B)). (118) 



They are highly constrained by the SU(N) x SU(2) symmetry. For example, 

(Bi) = (119) 

using SU(2) invariance. This might seem strange, because we were using an expansion 
Bi = Xi+Ai with Ai being a "small" fluctuation. However, this is not a contradiction: we 
do not assume any spontaneous symmetry breaking (or gauge fixing), and the solution 
Bi = Aj is only one possible gauge choice. To get nontrivial results, we should of course 
consider quantities which are gauge invariant or contain gauge invariant information. 
Consider for example 

(Fi^Fjiy)) = c 5ij g ab v a ® v h (120) 

using the SU (N) x SU(2) symmetry, where v a denotes the SU(N) Gell-Mann matrices. 
Here x and y stand for the first respectively second tensor slot, interpreted as functions 
on Sjj. Since g a b v a ® v b is the reproducing kernel, it should be interpreted as delta- 
function 5(x,y) on the sphere, and we can write 

(F i (x)F j (y)) = c6 ij 6(x,y). (121) 
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This makes sense: there are no propagating modes, therefore there is no correlation 
between fields at different points. The normalization can be calculated e.g. for coinciding 
points, J {Fi(x)F l (x)) which is certainly nonzero but finite. Similarly, it follows that 

{Fi{ x )Fj{y)Fk{z)) = e ijk f abc u a ®v b ®v c (122) 

where f a b c is the structure constant of SU(N), since the r.h.s. is the only invariant 
tensor (d a bc cannot occur since the lhs is totally symmetric). Hence it is proportional to 
the function 

f(x,y,z) :=f ahc v a ®v h ®v c (123) 

on Sff, which in the classical limit is the unique function f(x, y, z) on the sphere which 
is invariant under APD's and vanishes for coinciding points. 

In general, one would expect that all correlation functions of the field strength have 
a well-defined classical limit in this model, just like the partition function. In principle 
it should be possible to calculate them explicitly in terms of integrals over eigenvalues, 
and we expect no problem related to UV/IR mixing or renormalization. It would also 
be interesting to know whether it is possible to relate them to integrable models. These 
issues are left for further investigations. 



7 The deformed fuzzy sphere revisited 



Remarkably, we can repeat the same construction with a slightly different potential, and 
obtain a gauge theory on the g-deformed fuzzy sphere [21]. Consider 



= \ + \f a\ q) (124) 

where A^ , i = 1,2,3 are N x N hermitian matrices and er^ are the g-deformed sigma 
matrices, both of which are defined to be the Clebsches of the corresponding irreps 
of U q (su(2)). To simplify the notation, we shall sometimes omit the label (q), which 
is understood throughout this section. We also need the g-deformed invariant tensors 
e % l = (^(q))^ and g^ = g\j which can be found e.g. in The generators satisfy 

e«X i X j = i\ k , XiXjg ij = [N ~ + 1]q (125) 

l Z \ q N 



where [n] q = q _ q _ 1 , and 

= i[2} q 2a\ a'a 1 ,,,, = [3} q (126) 

which implies 

oV' = %e l la k + g ij . (127) 
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It follows that 

(A %)) - ~{\ cr (q)J^ m • \ LZH > 



This means that has eigenvalues ±, / — )M JV+1 ^ + 7 with multiplicities N ± 1. It 
is therefore a minimum of the matrix potential 

^--M 02 -' 1 ^^ 1 ^- (129) 

Expanding now a general matrix C as 

C=^+Bio{ q) , Bi = Xi + Ai (130) 
imposing again the constraint C = |, one obtains 

p2 _ ( [JV _ iy„ + llt + 1 } = ( w _ [N _ ^ + llg ) + a% (i3i) 

where 

= Bi + ie^BkBt, F kl := z[A fc , ^] - i[A,, A] + i[A k , A] + s klm A m (132) 



is indeed the appropriate g-field strength as used in [23]. Hence we naturally recover 
gauge models on the g-deformed fuzzy sphere. However it is not clear how to define the 
trace: If we take the classical trace 

S = TrV q (C), (133) 
a strange term of the form Tr ((BiB^F-^) appears, because cr^ = f g ^-1 J ^ as a 



H°ia)J M ) = 2g&+«q-q- 1 )el (134) 



non-vanishing trace: 

Hence it seems more natural to take the quantum trace over the a 1 space, which has the 
property that tr g (cr l , q \) = 0, so that t r q{ a \ q ) a \ q )) = [2] 9 g^y Then taking the classical 
trace over the N x N matrices would lead to the action 

2 / N 2 — 1 \ 
S = -^Tr^BiB* — ) 2 + (B t , + ie m &B k )(B i + ie irs B r B s )y (135) 

where g-deformed tensors are understood. This is again invariant under SU (N) (hence 
solvable), but the "physical" rotations defined similar as in Section El are violated. On 
the other hand, taking the quantum trace over the full 2N x 2N matrices would break 
the gauge invariance but make the model formally U q (su(2)) invariant as in [21]. This 
may shed some light on the issues raised in the application of g-deformed gauge theories 
on D-branes, see |37j . 
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8 Discussion and outlook 



We presented a new formulation of pure gauge theory on the fuzzy sphere, which allows 
to carry out the path integral explicitly. The partition function for U(n) Yang-Mills 
theory on Sfj is calculated in the large N limit, and the known result [3 EH for the 
classical sphere is recovered for large N. 

There are several messages that should be stressed. First and foremost, gauge theo- 
ries on noncommutative spaces are accessible to new methods and tools which could not 
be applied on commutative space. Of course, Yang-Mills theory on S 2 is a rather simple 
field theory with no propagating degrees of freedom; however it does have nontrivial 
monopole sectors, and its quantization is not entirely trivial. If the methods presented 
here can be generalized, noncommutative gauge theory may become a useful alternative 
to lattice gauge theory, even from an analytical point of view. Of course it should also 
be interesting from a numerical point of view. 

Another important message is that the "classical limit" of pure gauge theory on 
the fuzzy sphere is smooth, at least for the partition function. This is not obvious in 
view of the UV/IR mixing effects in noncommutative field theories. It would be very 
interesting to know if this generalizes to higher dimensions. An interesting such space 
is fuzzy CP 2 [SH], which is the subject of current investigations. 

Furthermore, this new formulation of gauge theory on the fuzzy sphere is arguably 
simpler than on the classical sphere. It is defined by a potential for a hermitian matrix 
plus a constraint. This leads not only to the correct kinetic terms, but also all the 
monopole sectors arise "automatically" in a very simple way, with the correct classical 
limit. In particular, there is no need to introduce nontrivial fiber bundles, connections 
and other mathematical structures in our approach. What is missing so far is the 
inclusion of fermions in this formalism; this will be discussed elsewhere. 

There are several other aspects which require further work. First, one should be able 
to calculate the correlation functions or other suitable observables for pure gauge theory. 
Perhaps there are some connections with integrable models. One may also try to extend 
this approach to other gauge groups. Of course it would be very desirable to simplify 
the calculation in Appendix B, and to systematically calculate the corrections for finite 
N. Furthermore, there exist other interesting gauge models with similar action but 
without the term {BiB % — n2 ^ 1 ) 2 , as discussed in |18j . They do not fit very well into the 
formalism presented here. However then the "radial" field ip (jUj) becomes dynamical, 
and these models describe branes on SU(2) near the origin rather than a Yang-Mills 
gauge theory on a sphere. All these questions certainly deserve further study. 
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Appendix A. Some useful formulae 

The irreducible N- dimensional representation of the su{2) algebra Aj (JBJ) is given by 

(A 3 )h = S k i N + \ — ~> ( 136 ) 
(K)m = S k+1 ^(N-k)k, (137) 

where k, I = 1, N and A± = Ai ± 1X2- 

Furthermore, recall that 

<rV S ij ■ ie ijk cr k . (138) 
Together with (JHJ), this implies the following crucial property of the matrix Aj<7 1 : 

N 2 - 1 

(A i0 r') 2 = — ^-(V), (139) 

which means that the eigenvalues (in Mat(2N)) of the matrix Ajcr* are ~ ^ . To get 
the multiplicities, we note that Ajcr* is an intertwiner of (2) ® (JV) = (N — 1) © (N + 1) 
(i.e. it is invariant under £77(2)), hence the multiplicities are N + 1 resp. N — 1. 



Appendix B: Evaluation of the partition function. 

Consider the expression 

^? 

in (|1()5|) . At first sight it may appear ill-defined, because the denominator is singular 
due to the form (jlOOj) of J. However, the fraction is analytic in J, because all poles are 
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canceled by zeros in the determinant (it must be so, because the lhs of (ll()5|) is clearly 
analytic). To see this explicitly, assume the following "regularization" 



for some infinitesimal constant e, and consider A(j7 7 ) in more detail. By treating the 
contributions from the 2 blocks as above separately, it is easy to see that 

A{Jj) = A 4 (A;,)e M . (142) 
where ki are the eigenvalues of K. It is useful to order the A as A = (A^ + \ A^ 

N 
2 



A[ \ A| \ An_ ) where A 4 - ~ ±4- We then have to evaluate the determinant 



det(e^)=det^ iA _ ft+e ; ^ J, (143) 

which clearly contains a factor e M due to the degeneracy in the {ki}. To proceed, we 
expand it by choosing M + n resp. M — n columns in the upper resp. lower block of 
(fTIHjl as follows 

det(e iA ^)= (-l) ff(U+} ' U " }) det(e iA ^ + )det(e iA ^") (144) 

{J+UJ-} 

where {J + }, {J~} C {ki + e, ki} are complementary subsets with |J7" + | = M + n, \J~~\ = 
M — n, and the sign is given by this choice of subsets. In the terms det (e lAi ^ ), the J + 
are assumed to be ordered as in (k± + e, ku + e,k±, /cm), so that only the choice of 
the subsets {J^ matters for the sign. In terms of the fluctuations 

N 

V] = (Aj + y), J = l,2,...,M-n, 
N 

H+ = (A+-y), j = 1,2, ...,M + n (145) 

this becomes 

det(e iA ^') = (-iy {{J+UJ ~ }) det (e^ + ) det (e^~) e^ J+ ~^ J '\ 

{J+UJ-} 

(146) 

Now the rapidly oscillating terms have been isolated in the last exponential, and it turns 
out that the correct expansion is in the number of these rapidly oscillating variables ki 
in e l ~ J . This depends on the choice of {J + } (which of course fixes {J~}): Let 
be the set of fcj's which occur twice in {J + } (for e = 0). Because \ = M + n, there 
are at least n such k + 's. Assume that |{k + }| = n + d: then there must be in addition d 
elements k~ among the {ki} which occur twice in {J~}. The last term in ()146|) is then 
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e iN(J2 « ) ( which means that there are n + 2g? rapidly oscillating variables among 
the ki. We can therefore expect that d = will give the dominating contribution, and 
concentrate on this case first. Then {J + } = {J~} U 2{/t + } for e = 0, and 

det(e iA ^)d=o = J2 (-iy aj+UJ ~ }) det (e^ + ) det (e^ r ) e iJV ^ *+). 

(147) 

We will omit the superscript of k from now on. Consider {J + } in more detail, for 
fixed {k}. It has the form {k} U {k + e} U {JT^ ± e}, and each different sign choice 
in {J~ ± e} gives a different contribution to the sum which must be added. Different 
choices are related by an exchange of the elements k^ ki + e ^ {k} in . This amounts 
to an exchange of the corresponding columns in (|143j) . hence the two contributions 
come with a relative — sign. Among these, there is one "ordered" choice (j7" + )o = 
(ki + e, km + e,Ki, K n ), (<7~)o = (k'i), which depends only on {k} and serves to fix 
the sign (— 1) CT ( K ) := (— iy((-7 + )o,(-?~)o) _ The following notation is useful: 

Tlf(k h k M ) = f(ki, .., h + e, k M ) (148) 

so that e.g. 

M 

(j + h = n T i(^ -> «i, , «n) = n ^ j+ ( i49 ) 

where J + = (fci, ku-, Ki, Kn), ^~ = k' M _ n ) from now on. Summing over all 

these permutations for each set {k} leads to 

det(e* A ^)d=o = ^(-l^V^*) ]J ((T l € )+ - det {e^ J t) det (e^ 7 )- 

{k} 

(150) 

Here (Tf)± indicates that the Tf operator acts only on det (e lfli j j ). The sign (— 1)°"( K ) 
now depends only on the choice of {k} C {ki}. Using 

(77)+ - m)- = mu - 1) - m)- - 1) - c((^)+ - W-) (i5i) 

for e — > 0, we get 

det(e iA< J >) d=0 = e M " n J^-l)' 00 e W(E k) \{ ((^)+ - det (e^ + J /) det (e^~ Jr ). 

{«} {fc'} 

We have hence recovered most of the e factors; the missing e n is contained in det [e % ^ 1 J j ) . 
To proceed, we can now integrate over ^ in fll()4j1 . noting that 

A (A) = N x A(/i + )A(/0 exp{n(EfT - E^ + )) (153) 
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up to corrections of order 1/N. In fact we can even neglect the last exponential compared 
to the leading terms exp(— ^ ^ (/x ± ) 2 ) in the action, expanded as in (JTBJ). Furthermore, 

j d^A(^) det(e^) eX p(-^ $>+) 2 ) = c A( J+) exp(-£^ J+) 2 ) (154) 

for some constant c by antisymmetry in J + , and similarly for jjT . Noting also that 

A(J+)A(J-) = ((-l) ,T(K) e n A 4 (/€)A(A; / ) ]J(k' - k) 2 ) + (a(A/)) (155) 
we get using (|142j) 

= / dh ^ y, ^ ex p ( iiV E «) - |^E * 2 )) 

n (w+ - W-) ( W - -) 2 ^') «p(-£^ E k ' 2 ))i^ k ') -p(-^ E fc 

always ignoring numerical constants. Notice that all signs have disappeared. Since 

(9 fcl -^ 1 )(/(fc 1 )/(/ 1 ))| fcl=il =0, (156) 



only the term Yl(k' ~ K ) 2 = (2? y -«■ ) IK^' — K ) 2 survives the derivatives w.r.t. fc^, 
and the term (£\ k ,^ K _ ) can be moved outside of the remaining derivatives 7 . Repeating 



1 yi^i t ' 

this, we find 



J d h e-^ £ A 2 (-)^ EK) exp(-^E^) (EII^) ( 157 ) 



since A 2 (£;) = A 2 (k)A 2 (A;') H(k' - k) 2 . Here 

iVi -j v, 



1 ^ 1 iV„ t 

1:11,.. ., 1:11,.. ., - ii,. 



is the sum over all possible splittings of {k'} into n partitions {k'}^, {k'}N n of size 
Ni + ... + N n = M — n, and flf* ^) = (ll {fc ' }jVi pz^r)- That is, each k[ is "linked" 
to one Kj by a factor of the form k , * K . , and there are Nj such fc' linked to Kj. For fixed 
(iVi, ...,N n ) each of these terms gives the same contribution, therefore we simply get a 

7 If we were more careful to include the factor exp(n(£/^~ — there would be additional 

contributions which are suppressed by 
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multiplicity factor jj^ Furthermore each different choice of {k} C {A;} gives the 
same integral since the form is identical. Therefore 

e ^fi(F^)---ri(F^)-(-^E* 2 )-(^ 9 ) 

A r i+...+A r „=A/-n n 

We have now taken just one fixed subset {k} C {k}, since all choices give the same 
result. 

The integral over the k[ has now the same structure for all i and can be carried out. 
Notice that there are poles in which seems a bit surprising because the original 

expression ^^x^nk is perfectly regular. However, recall that this corresponds to the 
sum over all possible terms in (|14fij) . in particular over all choices of {k} C {k}. In 
(jl59j) these different contributions are included in the integral over all values dki, 
contributing with the same \ki — Kj\ but opposite sign. Therefore the cancellations 
in the sum over {k,} are now reflected in the cancellations in the contributions to the 
integral from both sides of the poles with the same \ki — Kj\. We should therefore use the 
principal value of the integral -f dki in (|159j) . which is perfectly regular and well-defined 
because there are only simple poles. To put it differently, we can restrict the range of 
integration to the space of fcj with \ki — kj\ < e', say; this must give the correct result 
for e' — > 0. But this is just the definition of the principal value of the integral. In 
fact, notice that if there were no poles in the above formula, each integral J dk[ would 
produce an exponential factor of order e~ N ^ g2 , and Z would vanish for large N. The 
contributions from the poles will give the correct, finite result. 



Using the identity 



du 1 -f(u)= 1 -(l du f{u) - f{ - Uy ) (160) 
u 2 \./ u 



we get 



2 

/\2\ 



J- dk' ^— eM-ik' - ^(k 



- . - — . - / • / • -u^k sin(it) • sinh(^/cw)\ g 2 2 

e ik 2N k i du[ - ie u n k — ±-L - e lu — — - exp( — - — u j 



u u ) 2N 

Now - sin(u) ~ tt5(u), and ^ smh(j^Ku) ~ under the integral, hence 



1 n 2 2 

dk' exp(-ik' - 4r- T (k') 2 ) -> -m e - iK -^ K2 for large N. (161) 

k — k 2N 
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We therefore obtain 

,2 



Jd=0 — 

Ni+...+N n =M-n 



fAT n =M-n 1 "" n ' ^ 

(162) 

In terms of the integers rrii = N — (JVj + 1) which satisfy 

m,j = nN — M =: m, (163) 

the combinatorial factor is 

(M-n)\ ^ (M-rQ! 

JVkLJVJ ~ ((£ - l)!)n ^ ' 

as long as the rrii are small (which will be justified below), up to corrections of order ^. 
Hence we can drop this factor, and obtain 

z d=0 = Yl J d « A '( K ) ^ ex P(-y E ft 1 ~ jf))- (165) 



mi+...+in„=m 



This is now a perfectly nice integral. The ^ in the exponential can be neglected, and 
(jTUfiJ) follows. 

Several remarks are in order. 



• If one would similarly treat the case of more oscillating factors d > in ([146)1 . we 
would get similar formulas with k replaced by n + 2d variables k + , k~ , which come 
with oscillating terms e lNKi . The remaining analysis would be similar, with only 
N — n — 2d variables k' which contribute an integral as in (jlfilj) . This leads to 
an expression similar as in ([162)1 . however there are now not enough phase factors 
e~ tKi to cancel the rapid oscillations in e lN ^ Ki . Therefore this integral will have 
additional rapidly oscillating terms, which lead to an exponential suppression. 
Therefore the contributions Zd>o can be neglected for large N. 

• The approximation replacing the combinatorial factor in ([164)1 by a constant is 
justified, since only small \rrii\ contribute to the final result ([107)1 . 

• one can understand the above calculations, in particular the dominant contribution 
from the poles in ([159J1 intuitively as follows: Each contribution from the poles 
with (N\, N n ) corresponds to a "clustering" ({k'}^ ~ Ki), ({A;'}jv„ ~ kJ, 
which gives the dominant contribution to the integral over the ki. These clusters 
correspond precisely to the saddle-points discussed in Section which are also 
clusters of Ni + 1 eigenvalues of C. Hence the leading contribution comes form 
the saddle-points. 
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• This calculation could be generalized to the potential (f78j) . ordering the eigen- 
values as A = (A[ + \Ai + \...,A^ + \Af\...,A[~\A i i\...,AL ) ) where A, (±) « ±f 
and Af^ « |. This would allow to explicitly calculate the relative weights of the 
different topological sectors in (|108|) . which have been put by hand here. 

• with some effort, it should also be possible to compute the leading correction 
terms in i The relevant approximations are those in (|164|) and (|16ip . which can 
certainly be improved. The other approximations (taking only d — 0, ignoring the 
exponential in in ()153j) an the higher terms in (fTKjl ) apparently give corrections 
which are exponentially suppressed. An exact calculation is certainly desirable, 
but would require more sophisticated tools. 

• needless to say, it would be nice to simplify this calculation. 
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